Square lattice site percolation at increasing ranges of neighbor interactions 
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We report site percolation thresholds for square lattice with neighbor interactions at various 
increasing ranges. Using Monte Carlo techniques we found that nearest neighbors (N 2 ), next nearest 
neighbors (N 3 ), next next nearest neighbors (N 4 ) and fifth nearest neighbors (N 6 ) yield the same 
p c = 0.592 • • • . At odds, fourth nearest neighbors (N 5 ) give p c = 0.298 ■ ■ ■ . These results are given 
an explanation in terms of symmetry arguments. We then consider combinations of various ranges 
of interactions with (N 2 +N 3 ), (N 2 +N 4 ), (N 2 +N 3 +N 4 ) and (N 2 +N 5 ). The calculated associated 
thresholds are respectively p c = 0.407 • ■ ■ , 0.337 • • • , 0.288 ■ • • , 0.234 • ■ ■ . The existing Galam- 
Mauger universal formula for percolation thresholds does not reproduce the data showing dimension 
and coordination number are not sufficient to build a universal law which extends to complex lattices. 

PACS numbers: 07.05.Tp, 64.60.Ak 
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I. INTRODUCTION 

Calculating percolation thresholds has been an ongo- 
ing challenge for decades QlSlilELIHEl- While very few 
lattices allow an exact analytical calculation, large scale 
simulations have been very valuable to determine a large 
spectrum of them for both Bravais' 0, 0, and disor- 
dered |3 lattices. The drastic increase in computer ca- 
pacities has recently permitted the calculation of thresh- 
olds at rather high dimensions up to d = 13 for the hy- 
percube [4J, [lfj . In parallel, not much work has been de- 
voted to regular lattices with neighbor interactions which 
are not nearest neighbors (N 2 , von Neumann's neighbor- 
hood). Some scarce results are available for simultane- 
ous nearest and next nearest neig hbors (N 2 +N 3 , Moore's 
neie hborhood) HE). 

In this paper we report for the first time a systematic 
calculation of site percolation thresholds for the square 
lattice with neighbor interactions at successive increasing 
range. We consider the series of nearest neighbors (N 2 ), 
next nearest neighbors (N 3 ), next next nearest neigh- 
bors (N 4 ), fourth nearest neighbors (N 5 ) and fifth nearest 
neighbors (N 6 ). It should be stressed that for each one of 
the considered distance of interaction, all others are not 
active. For instance in the case of next nearest neighbors 
(N 3 ), the nearest neighbors (N 2 ) sites are not connected, 
only the N 3 are. This principle applies to all our calcu- 
lations. We found that the threshold is the same for all 
of them with p c = p c (N 2 ) except at N 5 . An explanation 
in terms of symmetry is provided. 

We then consider combinations of various ranges of 
interactions with (N 2 +N 3 ), (N 2 +N 4 ), (N 2 +N 5 ) and 
(N 2 +N 3 +N 4 ). In these cases we have simultaneous range 
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of interactions but they are necessary compact. For in- 
stance for (N 2 +N 4 ) all nearest neighbors sites are con- 
nected as well as all next next nearest neighbor ones but 
next nearest neighbors are not interacting. 

Comparing our numerical estimates with the predic- 
tions from the Galam-Mauger (CM) universal formula 
for percolation thresholds |l2(, we found significant dis- 
crepancies. It strengthens the earlier claim that only di- 
mension and coordination number could not be sufficient 
to build a universal law which extends to complex lattices 

01a. 



II. CALCULATIONS 

There exist several computational techniques which 
allow to perform calculations of percolation thresholds 
HH [IB EH 113 Here we are using the Hoshen-Kopelman 
algorithm (HKA) [l7|. Once the lattice is given with the 
occupied sites, it allows to recognize which sites belong 
to which clusters. With HKA one can assign to each 
occupied site a label and sites in the same cluster have 
the same labels. Different labels are assigned to differ- 
ent clusters. The HKA is particularly efficient when we 
check if the site at distance I from the first line — often 
fully occupied — is still connected to that line through 
the sites at the distances smaller than I. The algorithm 
requires storing only single line of sites and goes through 
the lattice only once. In such a case HKA becomes ex- 
tremely efficient as it saves memory and time Q. How- 
ever, when links between sites at distances larger than £ 
from a top border are desired, the whole lattice must be 
stored [j, EU |2(| . With the HKA on a square lattice 
when we assign the labels for the investigated site (black 
sites in Fig. [TJ, we need to check already labeled and 
occupied sites in its neighborhood (slashed sites in Fig. 

The possible links to remaining sites in the neighbor- 
hood (backslashed sites in Fig. ^) may be checked later, 
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FIG. 1: Various site neighborhoods on the square lattice: (a) 
N 2 — von Neumann's neighborhood, (b) N 3 , (c) N 4 , (d) N 5 , 
(e) N 6 , and the examples of their combinations: (f) N 2 +N 3 
- Moore's neig hborhood, (g) N 2 +N 4 , (h) N 2 +N 5 and (i) 
N 2 +N 3 +N 4 . 



basing on the neighborhood's point symmetry. 

The percolation thresholds values p c are evaluated 
from the crossing point of three curves showing depen- 
dences of the percolation probability P on the site occu- 
pation probability p for lattices of linear sizes L = 100, 
500 and 1000. The results are averaged over N TUa = 10 3 
and 10 4 for L — 1000 and 100, respectively. With enlarg- 
ing the lattice size L the curve P(p) becomes stepper and 
stepper and tends to Hevisade's function 0(— p c ) when 
L — > oo, as expected. 



III. RESULTS 

We present our results in Tab. [I] The percolation 
thresholds p c for the square lattice are computed with 
HKA for a series of neighborhoods. First only one type of 
neighbor interactions is considered at a time, increasing 
repeatedly the range with N 2 , N 3 , N 4 , N 5 and N 6 . It 
turns out that the threshold p c = 0.592 • • ■ is the same 
for all of them except at N 5 where p c = 0.298 

Indeed, all lattices with neighborhoods shown in Figs. 
Bib), (c) and (e) may be mapped into a N 2 situation as 
in Fig. [Ha). The only difference is a larger and larger 
lattice constant. To implement the mapping, we take a 
square lattice and build on it the lattice from only N 3 
interactions. Two independent interpenetrated squares 
sublattices appear. Therefore the percolation of N 3 is 
split onto two parallel N 2 problems on each one of these 
two square sublattices. Accordingly the p c on each one is 
the p c of N 2 . Moreover, as the site must be distributed 
homogeneously on the initial lattice, we will have the 



TABLE I: The percolation threshold p c for various neighbor- 
hoods on square lattice and sites coordination number z and 
the theoretical values p^ M . 



neighborhood 



N 2 


4 


0.592- •• 


0.5984 ■ ■ ■ 


N 3 


4 


0.592 ■ ■ ■ 


0.5984 ■ ■ ■ 


N 4 


4 


0.592- • • 


0.5984 ■ ■ ■ 


N 5 


8 


0.298- •• 


0.4411 ■ ■ ■ 


N 6 


4 


0.592- •• 


0.5984- •• 


N 2 +N 3 


8 


0.407 ■ ■ ■ 


0.4411 ■ ■ ■ 


N 2 +N 4 


8 


0.337- • • 


0.4411 ■ ■ ■ 


N 2 +N 5 


12 


0.234- •• 


0.3748 • • • 


N 2 +N 3 +N 4 


12 


0.288- •• 


0.3748 ■ • ■ 



(a) 



(b) 




FIG. 2: The lattices with (a) N 3 and (b) N 4 neighborhoods 
may be mapped to (a) two or (b) four parallel N 2 situations 
but with (a) s/2 and (b) two times larger lattice constants. 



same density of occupied sites on each one of the sub- 
lattice making both percolation to occur simultaneously 
at the same p c (see Fig. Efa))- Such a scheme can be 
repeated with N 4 (Fig. □») and N 6 but not with N 5 . 
As shown in Fig. ^d) the N 5 lattice has eight neighbors 
while N 2 , N 3 , N 4 and N 6 have four. 

These symmetry properties may become instrumental 
in underlining interesting physical properties associated 
to some exotic materials. In particular if one is able to 
discriminate between the two interpenetrated lattices, it 
may open a way to reach percolation at a much lower 
critical density, down to halt the value of p c . But such a 
search is out the scope of the present work. 

We also consider several combinations of various ranges 
of interactions. First an increasing compact neigh- 
borhood with (N 2 +N 3 ) and (N 2 +N 3 +N 4 ). The cal- 
culated threshold numerical estimates are respectively 
p c = 0.407 • • • and 0.288 • • • (Tab. CJ). Then more com- 
plex ones with (N 2 +N 4 ) and (N 2 +N 5 ) for which we ob- 
tained p c = 0.337 • ■ • and 0.234 ■ ■ • (Tab. [J). The fact 
that p c of (N 2 +N 5 ) is smaller than p c of (N 2 +N 3 +N 4 ) 
is consistent with N 5 z — 8 instead of z — 4 for all the 
others. 

The obtained percolation threshold values p c (N 2 ) = 
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0.592 • • ■ and p c (N 2 + N 3 ) = 1 - p c (N 2 ) = 0.407 • ■ • are 
consistent with the values reported in Refs. 0, 0, 0. 
We have also revised the value of p c (N 2 +N 3 + N 4 ) which 
has been studied in an earlier paper |22^ and it was put 
at 0.292 to compare to our value of 0.288 

IV. DISCUSSION 

At this stage it is interesting to check the validity GM 
universal formula for percolation thresholds [12J, |23| in 
the case of these complex neighbor interactions. Com- 
paring our numerical estimates with its predictions as 
shown in Tab. [I] we found a good agreement for N 2 , N 3 , 
N 4 and N 6 (A = 0.006) but not for N 5 (A = 0.123). It is 
also fair for (N 2 +N 3 ) with A = 0.034 but not for all oth- 
ers combinations. The significant discrepancies occur for 
complex and non compact neighborhood. It strengthens 
the earlier claim that only dimension and coordination 
number could not be sufficient to build a universal law 
which extends to complex lattices |l3[ . 

Indeed above failures could be anticipated due to the 
fact that several lattices have both identical z and d tough 
they exhibit different thresholds as seen from Tab. [I] In 
particular N 5 , (N 2 +N 3 ), (N 2 +N 4 ) have z = 8 and d = 2 
while all p c are different. The same occurs for (N 2 +N 5 ) 
and (N 2 +N 3 +N 4 ) with z = 12 and d = 2. 

Similar situation occurs for Tc in the Ising model 
where even with the same number of interacting spins in 
the neighborhood and the same dimensionality we have 
different Tq |24j. On the other hand Bragg- Williams 
approximation I25j predicts Tq to be unique function of 



coordination number z, i.e. ksTc = zJ 26]. The GM 
universal formula which also extends to Tc includes a 
dependence on both d and z |27| . 

To conclude we have report for the first time numerical 
estimates for site percolation thresholds for the square 
lattice with N 3 , N 4 , N 5 , N 6 , (N 2 +N 4 ) and (N 2 +N 5 ) 
interactions. Our new estimates may prove useful in 
the search for a robust universal formula for percolation 
thresholds which would apply to complex lattices. In par- 
ticular on how to extend the GM law by including some 
additional topological ingredient besides coordination z 
and dimension d. 

These results may prove useful to some of the large 
spectrum of physical and interdisciplinary topics where 
the percolation theory may be app lied like forest fires 
spreading |2Ct l28l . immunology 12 9|. liquid migration in 
porous media |3(j, econophysics [31|, and sociophysics 
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